Abstract 

We have studied the angular dependence of the irreversible mag- 
netization and its time relaxation in YBa2Cu307 single crystals with 
one or two families of columnar defects inclined with respect to the 
c-axis. At high magnetic fields, the magnetization shows the usual 
maximum centered at the mean tracks' orientation and an associated 
minimum in the normalized relaxation rate. In contrast, at low fields 
we observe an anomalous local minimum in the magnetization and 
a maximum in the relaxation rate. We present a model to explain 
this anomaly based on the slowing down of the creep processes arising 
from the increase of the vortex-vortex interactions as the applied field 
is tilted away from the mean tracks' direction. 

1 Introduction 

Pinning of flux lines by columnar defects (CD) in high temperature super- 
conductors has been of considerable interest in the last years. It is well 
known that these correlated defects yield a strong enhancement of flux trap- 
ping, in particular if the applied fleld H is aligned with the tracksfl], ||, ^, |]. 
When H is tilted away from the linear defect direction beyond a lock-in an- 
gle Ql, vortices form staircase structures with kinks connecting segments 
trapped in the columnar defects. The appearance of these kinks is ex- 
pected to reduce the critical current Jc and to produce a faster relaxation. 
Thus, the angular dependence of the persistent current density J, should 
show[^, ^ a peak at the CD's direction, as indeed observed in many 
casesi, i I, I, I, |10|, O, mil, 0, |T1 . 



However, Zhukov et al.[|^ have recently shown that the angular depen- 
dence of the irreversible magnetization of YBa2Cu307 (YBCO) crystals with 
CD along the c-axis exhibit a local minimum rather than a maximum for 
that fleld orientation. They found that this interesting and anti-intuitive 
behavior is related to geometrical effects; if the rotation axis is parallel to 
the shortest side of a rectangular sample the minimum is observed, but if 
the axis coincides with the longest side the usual maximum in the angular 
dependence is recovered. Their interpretation also involves a sharp increase 
in the critical current density parallel to the rotation plane as the fleld is 
tilted away from the CD. 
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Although the geometrical aspects of the anomalous behavior were con- 



vincingly demonstrated |T^, the origin of the increase of the current density 
as kinks proliferate is still very unclear. They speculate that it may be re- 
lated to the appearance of helicoidal instabilities in the kink structure, but 
certainly other explanations cannot be discarded based solely on their results. 
One important fact to be taken into account is that the persistent currents 
determined in magnetization studies of HTSC is usually much smaller than 
Jc, as it is strongly reduced by thermal relaxation. Thus, the observed fea- 
tures are more likely to be related to differences in the activation energy of 
the excitations that dominate the depinning process for different angles. 

In this paper we show that the anomalous dip is also visible in YBCO crys- 
tals with one or two families of aligned columnar defects inclined with respect 
to the c-axis. When only one family is present (all the defects are parallel) 
the local minimum is centered at the CD's direction. This result demon- 
strates that the anomaly is only due to vortex-tracks interactions, and the 
influence of crystal anisotropy or pinning by twin boundaries can be ignored. 
We also find that when two families of tracks are present (planar splay) only 
one minimum, centered at the mean defect's direction, is observed. Because 
in this case no kinks connecting pins of the same family are present in the 
angular range in between the two tracks' orientations, helicoidal instabilities 
are ruled out as a possible origin of the anomaly. To explore the nature 
of the thermal activation processes we performed time relaxation measure- 
ments as a function of angle. We show that the minimum in the irreversible 
magnetization is associated to a faster relaxation. We propose an alternative 
explanation of the anomalous angular dependence based on the reduction of 
the creep processes due to the increase of the vortex-vortex interactions as 
H is inclined with respect to the mean CD direction. 



2 Experimental 

We carried out magnetic studies of two YBCO single crystals grown from the 
self flux method and oxygenated following the procedures described in ref.l7. 
Both crystals were taken from the same batch, and display a Tc = 91. 6-?^ 
before irradiation. The two crystals have similar thickness t ~ 15/xm and 
approximately rectangular shape, with dimensions L x s ~ 0.67 x 0.22mm^ 
for sample A and L x s ~ 0.58 x 0.48mm^ for sample B, where L and s are 
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the long and short side respectively. 

Columnar defects were created by 309Mev Au'^'^^ ion irradiation at the 
TANDAR accelerator in Buenos Aires, Argentina. In both cases the defects 
were introduced |jl5| inclined with respect to the c-axis, with the irradiation 
plane (the plane formed by the c-axis and the irradiation direction) perpen- 
dicular to s. Sample A was irradiated at an angle 0d = 10° off the c-axis, 
and the dose was equivalent to a matching field = 3T. Sample B has two 
sets of columnar defects, one at Qdi = +5° and the other at Qd2 = +15°, 
each one with a matching field -B$i = -B$2 = 1.5T. In this way we obtain the 
same total dose (3T) and the same average angle for the columnar defects 
(10°) in both samples. 

DC magnetization measurements were made in a Quantum Design SQUID 
magnetometer with a 5T magnet. The magnetometer is equipped with two 
sets of detectors, which allows us to record both the longitudinal (Mi) and 
the transverse (Mf) components of the magnetization (parallel and perpen- 
dicular to H respectively). The samples can be rotated in situ around an 
axis perpendicular to H using a home-made rotating holder [|1^ . 

To perform the magnetic measurements the crystals were carefully aligned 
with the rotating axis normal to the irradiation plane, in such a way that the 
condition H || tracks could be achieved within ~ 1°. This configuration also 
satisfies the geometrical condition (rotation axis parallel to the short side) 
required to observe the minimum. 

It is known that the measurement of Mt in a QD magnetometer possesses 
some difficulties arising from the presence of a spurious signal due to the 
longitudinal component M; that is detected by the transverse pick-up coils. 
This occurs when the sample is slightly off-center with respect to the vertical 
axis of the coils, which is frequently the case. We have completely and 
satisfactorily solved this problem. The solution includes an initial alignment 
procedure and the external processing of the original SQUID output signal 
using software developed ad-hoc. All the details related with the hardware 



and software of the sample rotation system will be presented elsewhere [18 



We performed isothermal magnetization loops maintaining a fixed value 
of the angle 9 between the normal to the crystal (c-axis) and the apphed 
field direction, and recording both components Mi{H) and Mt{H). We use 
the widths of the hysteresis AMi{H) and AMt{H) to calculate the modulus 
Mi = ^JAM^ + AMf and direction of the irreversible magnetization vector 
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Mi. Loops were recorded up to if = 5T in all cases. As H is reduced 
from this maximum field, the non-equilibrium currents that generate the 
critical state profile start to reverse direction |T^. The formation of a fully 
developed critical state of the oposite sign occurs after a field decrease of 



the order of AH ~ 2H* ~ Jt, where H* is the self field 201. This situation 



is clearly identified as Mi(H) and Mt{H) reach the field-decreasing branch 
of the loop. We have carefully checked that all the AM data shown in 
this work correspond to the difference between two oposite fully developed 
critical states, thus Mj can be easily related to the persistent currents. After 
each loop is finished the sample was rotated, warmed up above Tc and then 
cooled down to the working temperature in zero field. In this way, the initial 
Meissner response was recorded for each angle. 

As the non-equilibrium currents in thin samples are strongly constrained 
to flow parallel to the sample surface. Mi points almost perpendicular to the 
surface 0, ^ |22|, ^ in a wide angular range of applied field < 6 < 0c. 
For both crystals the critical angle Gc ~ arctan(L/t) > 87°, and we indeed 
confirmed that Mi was normal to the sample surface within our 1° resolution 
for all the angles shown in the present work. (The angle G was determined 
independently using the Meissner slopes, as described in ref. 21). 



3 Results and Discussion 

Figure 1 shows the angular dependence of the modulus of the irreversible 
magnetization as a function of G for crystal A at two temperatures. The 
main feature of this figure is the evident asymmetry with respect to the c-axis, 
which is due to the uniaxial vortex pinning produced by the inclined columnar 
defects. The anomalous minimum is apparent at both temperatures. This 
dip is centered at the tracks' direction Qd = 10° (except at very low fields, 
as discussed below). At T = 35i^ the minimum is visible for all values of 
the applied field. Its depth first increases with H, reaches a maximum at 
if ~ IT and then progressively decreases. At T = 70i^, on the other hand, 
the dip is only observed at low fields, its depth monotonically decreasing with 
H until the behavior switches to the well-known peak at higher fields. At 
this temperature it becomes clear that the dip is " mounted" over the broader 
usual maximum centered at the tracks' direction. The angular width of the 
minimum decreases with both temperature and field increase. 
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Also visible in figure 1 is a shift of the dip from the tracks direction 
towards the c-axis. This shift occurs for the lowest fields at both temper- 
atures, although it is not shown at 70K for clarity. In a previous work we 
have shown that this effect is related to the misalignment between the inter- 
nal flux density B (which represents the vortex direction) and H, due to the 
anisotropyjp. From now on we will concentrate on the field regime where 
B II H. 

As a first step to investigate the origin of the anomalous minimum, we 
must determine the geometrical relation between Mi and the non-equilibrium 
currents flowing in our samples, as a function of 0. To that end we will use 



the extended Bean critical state model for in-plane anisotropic currents p4 
For our thin and approximately rectangular crystals we obtain 

M. = ^(l-— ^) (1) 
20 ^ 3L J2 

where Ji{Q) and ^2(0) are the current densities (constrained to flow in the 
plane of the sample) parallel to the long and short sides of the sample respec- 
tively, as sketched in figure 2(a). Eq. 1 is valid provided that J1/J2 < L/s. 

When H is parallel to the c-axis all the currents are perpendicular to H, 
thus the Lorentz force on the vortices is maximum for both current directions 
and then Ji(0 = 0) and J2(0 = 0) are equal to the critical currents Jd and 
Jc2 respectively. Note that, in contrast to the case analyzed by Zhukov et al., 
</ci(0 = 0) and Jc2(0 = 0) in our case are different, due to the inclination of 
the tracks with respect to the c-axis 1 2^ . 



We must now analyze how Mj is expected to behave when H is tilted from 
the c-axis. In this case ^2(0) remains perpendicular to H and thus ^2(0) = 
Jc2(0). On the other hand, only the component of Ji(0) perpendicular to H 
contributes to the Lorentz force p6|, and thus Ji(0) = Jci(0)/cos(0). Thus, 



if Jci and Jc2 were independent of 0, the Mj(0) would increase as grows 
from to 90°. However, the minimum originated in this effect is centered at 
the c-axis and not at the direction of the columnar defects. Moreover, the 
observed minimum in Mj is much sharper than l/cos(0). 

We then conclude that, also in our inclined defects' configuration, any 
explanation of the minimum must involve an increase of either Jci(0) or 
Jc2(0) as H deviates from the tracks' direction. As the second possibility 
has been ruled out by the results obtained when the crystals are rotated 
around the longer axis|jT6|, we will focus our analysis in the possible reasons 
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for the sharp increase of Jci(0)- 

The interpretation suggested by Zhukov et al. is based on the appearance 
of kinks connecting nearby tracks when the field is inchned with respect to 
them. When the sample is rotated around its shorter axis, those kinks are 
on the average perpendicular to J2 and have a component parallel to Ji. It 
was speculated that in the force free configuration associated to Ji kinks may 
develop helicoidal instabilities|2^, thus resulting in an increase of Ji. 

To check this possibility, we repeated the study on crystal B. The key 
difference in this case, as we will demonstrate below, is that for any field 
orientation in between the two families of tracks, kinks connecting tracks of 
the same orientation do not exist, and consequently helicoidal instabilities 
cannot develop. 

Figure 3 shows Mj(B) at T = 60K for crystal B. The anomalous mini- 
mum at low fields is also clearly visible in this case, switching to the usual 
maximum at high fields. Now both the minima and the maxima are cen- 
tered at the mean tracks' orientation, 6 = 10°. These Mj(0) curves cannot 
be satisfactorily adjusted by superposition of two minima (or maxima) cen- 
tered at Qdi and 0d2, indicating that the observed behavior results from 



the combined interaction of vortices with both families of tracks 15 



We now analyze the vortex structure in this crystal as a function of 9. 
For 5° < < 15° vortices may zigzag between tracks of different families. 
If two such tracks physically intersect, no kink is required to connect the 
two pinned vortex segments. If, on the contrary, the two tracks are not in 
the same plane, a kink connecting both pinned segments must exist. Let 
us consider a track of the family Qd2 = 15°. The number of tracks of the 
other family (0di = 5°) that approach to it within a distance D of its 
axis is equal to the number of such tracks that cross the rectangle of area 
A = 2L)t[tan(15°) - tan(5°)] ~ OMtD, as seen in figure 4(a). 

We can estimate such number as = AB^2/^o- This gives, for instance, 
an average of 23 "close approaches" within a distance D = 6nm, which ap- 
proximately corresponds to the diameter of the tracks. The average distance 
between such crosses is about S t/n ^ 700nm. The same estimate for 5 
can be obtained using the more elaborate analysis of Hebert et al.[^. 



The above analysis indicates that for field orientations G/ji < B < 
and at low fields, when vortex-vortex interactions are small, the energetically 
most convenient configuration for most of the vortices is to zigzag among 
tracks that intersect within their diameter {D < 6nm), thus not forming 



6 



kinks. Of course, considering the random distribution of tracks there is a 
probabihty that some kinks will exist, but the number of them will be neg- 
ligible as compared to the case of one single family of tracks at similar field 
and inclination. 

The situation in this angular region changes as H increases, because the 
vortex-vortex interactions tend to inhibit the transverse displacements re- 
quired to zigzag. For instance, for a field direction 6 = 10°, pinned segments 
of length 6 imply transverse displacements of ~ 6sin{5°) ~ 60nm. When the 
distance between vortices decreases to around this value (that corresponds 
to if ~ 0.5T), it becomes energetically convenient to form some kinks con- 
necting tracks of different families that do not intersect (i.e., separated by 
D > 6nm). This reduces S cx 1/D, and consequently the transverse ex- 
cursions. However, as long as D remains smaller than the average distance 
between tracks of the same family, d ~ 36nm in our sample, it will still be 
convenient to form kinks between tracks of different families (inter- family 
kinks) rather than within the same family (intra-family kinks). 

These inter-family kinks cannot develop helicoidal instabilities, for sev- 
eral reasons. First, their lengths are in the range of ~ lOnm, too short to 
entangle. Second, according to Indenbom et al. [27] these instabilities are 
only visible in extremely low pinning crystals, which is certainly not our case. 
Finally, the kinks have their main component perpendicular to the plane of 
irradiation, which in our geometry means parallel to the axis of rotation. 
Thus, the force free configuration required for the appearance of the insta- 
bilities may only be produced by J2, instead of Ji, contrary to the original 
argument. Consequently, the observation of the minimum in this angular 
range, clearly seen in figure 3, rules out the possibility that it is associated 
with the helicoidal instabilities. 

For 6 < 9di and > Od2 the nature of the vortex structure changes. 
It will now be energetically favorable for vortices to form staircases with 
segments pinned mainly in tracks of one family, connected by kinks of the 
same type and orientation as those formed in samples with all pins in a single 
direction. Therefore, to a first approximation we can ignore the second family 
of tracks, and we have a situation similar to that of sample A. However, it is 
clear from fig. 3 that no hint of a change in behavior is seen either at = 5° 
or = 15°. Thus, the minimum appears to be independent of the presence 
or absence of intra-family kinks. 

Once the helicoidal instabilities have been discarded as possible sources 
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of the anomalous minimum, we must search for an ahernative explanation. 
A fact that we have not considered up to now is that, due to the large influ- 
ence of thermal fluctuations on the vortex dynamics in HTSC, the persistent 
current density J determined through dc magnetization measurements in the 
typical time scale of SQUID magnetometers is much smaller than the "true" 
critical current density Jc- This suggests that the anomalous minimum may 
be related to the angular dependence of the time relaxation of J. 

To confirm this possibility we have measured the normalized time re- 
laxation rate of the irreversible magnetization, 5* = —dLn{Mi) / dLn{t) , for 
the splayed sample 5 as a function of 0. Measurements were performed at 
T = 60K for two values of field: H = 3T, where Mj(G) shows the usual 
maximum at the mean tracks' direction, and H = 0.5T, where Mj(G) ex- 
hibits the anomalous minimum. The curves S{Q) are presented in figure 5(a) 
and 5(b) respectively, together with the corresponding Mj(6) data, already 
shown in fig. 3. 

Before we discuss the data shown in fig. 5, we must analyze how the 
quantity 5* defined above relates to the normalized relaxation rates of both 
current densities fiowing through the crystal, 5*1 = —dLn{Ji) / dLn{t) and 
^2 = —dLn{J2)/dLn(t). Operating on eq. (1), we obtain 



5(6) = (1 - K{e))SiiQ) + K{e)S2{e) (2) 

Thus, the infiuence of the relaxation rate of Ji and J2 on the global rate 
is weighted by the angle dependent geometrical factor 

^(e) = . (3) 

6LJ2 — sJi 

The range of validity of Eq. 2 is the same as eq. 1, i.e. sJi < LJ2, or 
equivalently K < 1/2. The condition K = 1/2 corresponds to the change in 
the shape of the inverted roof, from that of fig. 2(a) to that of fig. 2(b), that 



has been identified|T^] with the maximum in Mj(B). On the other hand, 
K is minimum at B = 10°, where J1/J2 reaches its smallest value. If we 
estimate that Ji{Q = 10°) ~ 73(6 = 10°), then K{Q = 10°) ~ 1/8. Thus, 
in the angular range of the minimum, S mostly refiects the behavior of Si. 

The high field behavior shown in figure 5(a) is in agreement with the 
theoretical expectations 0: S decreases as the field orientation approaches 
the tracks' direction, due to the growth of the pinned fraction of the vortices 



8 



and consequent increase of the activation energy. The observation of a single 
minimum at the mean direction of the tracks is well described by the scenario 
discussed by Hwa et al.[^], according to which the forced entanglement of 
the vortices in the angular range 5° < O < 15° tends to inhibit the thermal 
relaxation. 

Figure 5(b) shows a different behavior. Here we observe a narrow and 
small peak mounted over a larger minimum, both centered at 10°. The main 
minimum (which is wider than aX H = 3T) corresponds once again to the 
increase of the pinned fraction. The central peak, on the other hand, is 
a new manifestation of the anomalous behavior. As we showed above, in 
this angular region S is basically a measure of Si. We conclude that the 
anomalous increase of Ji as H is tilted away from the mean tracks' direction 
is a consequence of the reduction of Si. It is important to note that the 
normalized relaxation rate is a very fundamental parameter of the vortex 
dynamics that characterizes the pinning and creep regimes, and is rather 
insensitive to the pinning details^]. On the contrary, the persistent current 
density is a more derived variable that depends on the time scale of the 
measurement. Thus, the basic concept is that the minimum in Ji(0) is a 
result of the maximum in S'i(G), and not the other way around. The goal 
now is to find the reason for this unexpected behavior of 5*1. 

It is well known that the increase of vortex-vortex interactions usually 
results in a decrease of the normalized relaxation rate, which manifests in 
a larger glassy exponent /i of the collective creep regime as compared to 
the single vortex creep. This stiffening of the vortex matter due to elastic 
interactions is very general, and rather independent of the pinning details, 
so it occurs both for correlated and random disorder. This suggests that 
the observed decrease in Si as we tilt the field away from the mean tracks 
direction may be due to the increase of the interactions. Strong support for 
this interpretation arises from two distinctive features of the M((B) data. 

First, we note in fig. 6(a) that, right at the mean tracks direction where 
the minimum occurs, Mj at low fields is independent of H. This is clearly 
seen in the inset of fig. 6(a), where Mj(8.4°) for crystal B at T = 60K 
is plotted as a function of H. The field independent Mj regime, which is 
characteristic of a system of non-interacting vortices, extends up to ~ IT, 
i.e., it roughly coincides with the field range where the minimum occurs. 

In the second place, if we tilt the field away from the tracks we observe 
that the increase of Mj(0) is steeper the higher H is. This means that in 
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the proximity of the minimum M{H) (at fixed G) grows with H, i.e. M{H) 
exhibits a fishtail shape. Fishtail loops (observed in many HTSC compounds) 
have been attributed to a variety of origins. In some cases BO] the increase of 
M with H has been shown to originate in the reduction of the relaxation rate 
with increasing H, which is a consequence of the increase of the vortex- vortex 
interactions. The increase of Mj with both the tilt angle AG = G — 10° and H 
suggests a common origin of both dependencies. This becomes apparent in fig. 
6(b), where the Mi(G, H) data of fig. 3 is replotted as a function of the field 
component perpendicular to the mean tracks direction, H± = Hsin{AQ). In 
the field range H < IT we observe that the various curves have the same 
curvature around the minimum. 

In summary, the scenario that emerges from the angle dependence of 
5*1 and the angle and field dependence of Mj is the following. In the field 
range of the anomalous minimum, and for H parallel to the mean tracks 
direction, vortex-vortex interactions are small. Those interactions increase 
with H_i, thus resulting in a reduction of 5*1 and the consequent increase of 
Mi measured at fixed time. 

The same features in Mi(G) are observed at T = 70K in crystal A: 
Mj(G = 10°) is field independent up to ~ IT, and the curves Mj(G) at 
different fields have the same curvature around the minimum when plotted 
as Mi{H±). 

The reason for the increase of the interactions with G at constant H, i.e., 
at constant average distance between vortices, is not obvious. In the case 
of zigzagging or staircase vortices, the distance between neighbors varies 
along the field direction. As the vortex-vortex repulsion is a highly nonlinear 
function of their separation, the strength of the interactions depends not 
only on their average distance but also on the amplitude of the transverse 
displacements. A detailed analysis of the interaction energy as a function 
of G thus requires the complete computation of all the 3D configurations 
involved. This is a very difficult problem, that we will not attempt to solve 
here. However, we will now present a simple estimate that shows how the 
interactions in the presence of splay defects and at low fields may increase as 
H is tilted away from the mean tracks' direction. 

As discussed above, for 5° < G < 15° and low fields, it is useful to consider 
that vortices zigzag within a planar grid of tracks as shown in fig. 4(b). The 
maximum displacement of the vortex perpendicular to the field direction, R, 
as a function of G, can be estimated as ~ 6sin{5°+ \ AG |). This relation 
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shows that i?(B) is minimum at the mean tracks' direction B = 10°. As a 
result, the minimum distance between two adjacent vortices decreases with 
A0, even when the average distance between them remains constant, thus 
producing an increase of the interaction. 

For sample A, with only one family of tracks, the transverse displacements 
must be calculated differently. Now the zigzag-vortices, without kinks, must 
be replaced by staircase-vortices with kinks connecting parallel tracks. The 
orientation of a kink (the angle between the kink and the tracks) depends 
on the pinning energy of the two adjacent tracks, thus the dispersion in the 
pinning energy of the columnar defects results in a dispersion of the kinks 
orientations [0. The larger the pinning energy, the closer to the ab-plane is 
the kink. In a previous study we have shown that, when exceeds the angle 
of a particular kink, such kink disappears and the vortex involved becomes 
trapped by stronger tracks connected by a longer kink, closer to the ab-plane. 
This process generates a progressive increase of both the average kink length 
and the deviation from B^) as AB grows, which again results in a decrease 
of the minimum distance between adjacent vortices, at constant H. 

We can use the analysis presented in the two previous paragraphs for 
the crystals B and A, with and without splay respectively, to compare the 
situation for B = 10° in both cases. At this angle, the average transverse 
displacements of the zigzagging vortices of crystal B is much larger than in 
crystal A, where vortices are expected to be locked. Thus, the interactions 
at the same H should be larger for crystal B. This is consistent with the 
observed values of the field required to switch from the anomalous minimum 
to the maximum in Mj(B). Although such field for each sample decreases 
with temperature, it is higher in sample A at 70 K (~ 3.5T, see fig. 1(b)) 
than in sample B at 60 K (~ 1.5T, see fig. 3). 

4 Conclusion 

The similarities in the behavior of the crystals with parallel and splayed 
tracks indicate that the physics involved in the anomalous minimum is rather 
independent of the details of the vortex configurations. The local maximum 
in the angular dependence of the normalized relaxation rate demonstrates 
that the minimum in Mj(B) is due to a stiffening of the vortex matter as H 
is tilted away from the mean tracks direction. We attribute such effect to 
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the increase of the vortex- vortex interactions arising from the enlargement of 
the transverse vortex displacements, but clearly further studies are required 
for a complete understanding of the phenomenon. 
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Figure 1: Irreversible magnetization Mj for crystal A (with a single family 
of tracks), as a function of the apphed field angle, 0, at several fields and at 
temperatures (a) T=35K and (b) T=70K. 

Figure 2: Sketches of the critical state profiles for different angular regimes. 
The pictures (a) and (b) shows the shape of roof pattern for ^ < j and 
^ > -J, respectively. The rotation axis is parallel to the shortest side of the 
sample, and perpendicular to the irradiation plane. 

Figure 3: Angular dependence of the irreversible magnetization Mj(if) at 
several fields for sample B (with splayed defects), at T=60K. 
Figure 4: (a) Sketch used to estimate the number of defects of family 1 
that approximate to one track of family 2 within a distance D, in a crystal 
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of thickness t. (b) Transverse excursions i? of a vortex in a planar grid of 
splayed columnar defects, function of angle. 

Figure 5: Angular dependence of the normalized relaxation rate S (open 
symbols) and irreversible magnetization M,i{H) (solid symbols) for sample B 
(with splayed defects), at T=60K and fields (a) H = 30kG and (b) H = 5kG. 
Figure 6: (a) Blow up of the data showed in fig.3 in the region of the minimum 
at low fields, (b) Irreversible magnetization Mj as a function of the field 
component normal to the tracks' direction, H±. The inset shows Mi as a 
function of the applied field H at the mean tracks direction. 
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